Coherent imaging of a pure phase object with classical incoherent light 
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By using the ghost imaging technique, we experimentally demonstrate the reconstruction of the 
diffraction pattern of a pure phase object by using the classical correlation of incoherent thermal 
light split on a beam splitter. The results once again underline that entanglement is not a necessary 
feature of ghost imaging. The light we use is spatially highly incoherent with respect to the object 
{~ 2/im speckle size) and is produced by a pseudo-thermal source relying on the principle of near-field 
scattering. We show that in these conditions no information on the phase object can be retrieved by 
only measuring the light that passed through it, neither in a direct measurement nor in a Hanbury 
Brown- Twiss (HBT) scheme. In general, we show a remarkable complementarity between ghost 
imaging and the HBT scheme when dealing with a phase object. 

PACS numbers: 42.50.Dv,42.50.Ar,42.30.Va 



I. INTRODUCTION 

Ghost imaging is a technique which allows to perform 
coherent imaging with incoherent Ught by exploiting the 
spatial correlation between two beams created by, e.g., 
parametric down conversion (PDC) E, 13 , B . B la l a B, 

ijijiMijiijiaB El m EE imif iflifii 

l23l l24 l25l l26l l27l l28l|7 Each of the correlated beams 
are sent through a distinct imaging system, traditionally 
called the test and the reference arm. In the test arm 
an object is placed and information about the object is 
recreated from the spatial correlation function between 
the test and reference arm. 

One of the most striking features of the ghost imaging 
technique is that since the two beams are spatially in- 
coherent no phase sensitive information about an object 
can be extracted by observing a single beam only. This 
means that the diffraction pattern of an object that sub- 
stantially alters the phase of the incoming light cannot 
be observed in any way in the test arm. Nonetheless, be- 
cause of the mutual spatial correlation between the two 
beams, we will show that the diffraction pattern may be 
reconstructed through the spatial correlation between the 
two beams. In other words, despite the beams being in- 
coherent, the coherence between them allows to perform 
coherent imaging: the scheme is therefore capable of do- 
ing coherent imaging with incoherent light. 

On the other hand, the diffraction pattern of an object 
that only alters the amplitude of the light, such as a 
Young's double slit, can be extracted from the intensity 
distribution of the object arm only, even when the object 
is illuminated with a spatially incoherent beam. One may 
do so by measuring the autocorrelation of the transmitted 
field as observed in the far zone of the object. This can 



'Corresponding author: lbache@coin.dtu.dkl Present address: 
COM»DTU, Technical University of Denmark, DK-2800 Lyngby, 
Denmark 



be conveniently done by using, the Hanbury Brown- Twiss 
(HBT) scheme |23- Thus, the ghost imaging technique 
is not required in this case, and one may stick to using a 
single beam only. 

It is therefore interesting to look beyond the case of 
an amplitude-only object. As an extreme case, in this 
work we want to observe the diffraction pattern of a pure 
phase object, an object that only alters phase information 

Emin. 

Initially the possibility of performing ghost imaging 
was ascribed to the presence of spatial entanglement be- 
tween the two arms |E Hi IM Hi- Lately this view has 
been challenged bymany groups lip. El Ellll IHIll 
ElESIMIlllilMllElllliiill (to cite but a few). 
Our group in particular has produced numerous publica- 
tions, showing both theoretically [13 and experimentally 
[l7l[l9ll27l | that basically all features offered by entangled 
ghost imaging can be mimicked by using a proper scheme 
that exploits classically correlated beams: the correlation 
between the beams is in this scheme created by dividing 
an incoherent pseudo-thermal speckle beam on a beam 
splitter (BS). The two outgoing beams are then still inco- 
herent on their own but since they are (classical) copies 
one of each other, they have a high mutual spatial co- 
herence. We showed that the only feature that cannot 
be mimicked by classical correlation is the 100% visibil- 
ity of the information, that can be in principle achieved 
only with entangled photons; however, in the classical 
case the visibility is still good enough to effecti vely re- 
construct the information in ElEl El millll 113. 
An important outcome of our analyses is that the entan- 
gled ghost imaging and our classical ghost imaging have 
the common feature of providing coherent imaging using 
incoherent light [l3|. Thus, both schemes should be able 
to reconstruct the diffraction pattern of any object, alter- 
ing amplitude as well as phase [l3|. We experimentally 
confirmed this prediction in the case of an amplitude ob- 
ject [13. 

An experiment of Abouraddy et al. ^ demonstrated 
the reconstruction of the ghost diffraction pattern of 
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a pure phase object using the entangled photon pairs 
produced by spontaneous parametric down-conversion, 
which represents an optical field that lacks second-order 
spatial coherence but is endowed with higher order spa- 
tial coherence. The introduction of Rcf. QJ claims that 
our ghost-diffraction experiment with classical thermal 
light reported in Ref. [I3| could be performed only be- 
cause effectively the thermal light was endowed with 
second-order spatial coherence. This might suggest that 
coherent imaging with spatially incoherent light is not 
possible with split thermal light. This interpretation of 
our experiment |17J | was not correct, because the light 
we used there was indeed spatially incoherent to a high 
degree. The viewpoint of Ref. jsl was challenged in 
Refs. [1^123, where not only we experimentally demon- 
strated again the ghost diffraction of an amplitude object 
with classical incoherent light, but also showed that the 
spatial incoherence of the light is a necessary ingredient 
to carry out the task. In this work we will finally demon- 
strate experimentally and theoretically that the claim of 
Ref. J3l is not correct also in the case of a pure phase 
object. 

We will show that the diffraction pattern of a com- 
monly used pure phase object, a transmission grating 
beam splitter (TGBS), can be reconstructed from the 
classical correlations between two highly spatially inco- 
herent beams generated by splitting a speckle beam on a 
BS. In order to render the light impinging on the object 
incoherent with respect to the object, we have to produce 
speckles of size on the order of 2.0 /im, and we achieve 
it through the so-called near-field scattering technique 
|30l l31| . Incidentally we remark that such a small speck- 
les size, which imply a spatial resolution on the order 
of 2.0 /xm, (see, e.g., Ref. [i3|) has no precedent to our 
knowledge in ghost imaging schemes, either with thermal 
or entangled beams. 

With such a degree of spatial incoherence, we will ver- 
ify that no information about the phase object diffraction 
pattern is present in the test arm, neither in the far zone 
intensity distribution nor in its autocorrelation, which is 
equivalent to a HBT type of measurement. The informa- 
tion gradually appears in the test arm as the degree of 
spatial coherence increases. The converse holds for the 
ghost diffraction scheme: the information on the phase 
object can be retrieved from the correlation between the 
test and the reference arm only when the light is spa- 
tially incoherent, and it disappears when increasing the 
coherence. From these results we will conclude that the 
claim of Ref. y\ is incorrect for the ghost imaging scheme, 
which indeed works as a coherent imaging scheme only 
because of incoherence. This claim could be possibly ap- 
plied to the HBT scheme, which in the case of a phase 
object works as a coherent imaging system only when the 
light is coherent. 

Thus, there exist remarkable differences between the 
ghost imaging and the HBT schemes, which will be clari- 
fied for the first time - to our knowledge - in the present 
work. 
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FIG. 1: Experimental setup for ghost diffraction of a pure 
phase object using classical pseudo-thermal light. See text 
for details. 



The paper is organized as follows. Sec. E] describes 
the experimental setup. Sec. lIIIl is devoted to theoretical 
results, with Sec. IIIlAl giving the basics of the theory of 
the ghost diffraction. Sec IIII B1 discussing the properties 
of the chosen object, and Sec. IIII Cl presenting some nu- 
merical results used as a guideline for the experimental 
implementation. In fact, the speckles needed are so small 
that we have to use a near-field scattering technique to 
produce them. This technique, along with its experimen- 
tal implementation is described in Sec. IIV Al Sec. IIVBI 
finally presents the experimental results of ghost diffrac- 
tion and of the HBT scheme with incoherent illumina- 
tion, while Sec. IIV (^ present results illustrating the tran- 
sition from incoherent to coherent illumination. Sec. Ivl 
concludes. 



II. THE SETUP 

The ghost diffraction setup used in our experiment is 
schematically shown in Fig. ^ A pseudo-thermal source 
generates a chaotic speckle beam, which enters a bal- 
anced 50/50 BS (the side of the glass cube is 12.5 mm). 
The pseudo-thermal source [Fig. H^a)] consists of a large 
coUimatcd laser beam (frequency doubled Nd:YAG laser, 
A — 0.532 /im, diameter D fn 10 mm) illuminating a 
slowly rotating ground glass followed by a square cell 
containing a colloidal turbid solution. The transversal 
size of the source is delimited by a pinhole of diameter 
Dph — 4.5 mm, which is placed at the exit face of the 
turbid cell. A detailed description of the thermal source 
together with the features of the speckles that it gener- 
ates will be given in Sec. IIV Al 

The BS divides the pseudo-thermal radiation into two 
"twin" speckle beams: the transmitted one is used for the 
test arm and is sent onto the object located right after 
the BS at a distance z = 18 mm from the source. The 



reflected beam is used for the reference arm; the mirror 
M deflects it towards the detector in a direction forming 
a smaU angle with the test arm. Note that, thanks to the 
double reflection of the reference arm, the test and ref- 
erence beams on the detector are not mutually reversed, 
but are each others replica. 

Both beams are collected by the central part of the 
lens F (focal length F = 50 mm) and their intensity 
distributions are detected by a charged coupled device 
(CCD) sensor placed in the focal plane of the lens. Due to 
the small angle between the two arms, the corresponding 
light spots fall onto different non-overlapping regions of 
the CCD sensor. 

The data are acquired with an exposure time « 0.1 ms 
much shorter than the speckle coherence time Tcoh ~ 
10 ms (see Sec. IIVA|) . allowing the recording of high 
contrast speckle patterns. The frames are grabbed at 
a rate of 5 Hz or smaller, so that each data acquisition 
corresponds to uncorrelated speckle patterns. 



III. THEORETICAL RESULTS 
A. General theory behind the setup 

In this Section we present the basic theory behind the 
ghost imaging setup of Fig. ^ 

The general theory of ghost imaging has been ex- 
plained in detail in Refs. [isL |27| . We summarize here 
the main points: 

1. The collection time of our measuring apparatus is 
much smaller than the time Tcoh over which the 
speckle fluctuate. Hence all the field operators are 
taken at equal times, and the time argument is 
omitted in the treatment. 

2. The speckle beam generated by the pseudo-thermal 
source is described by a thermal mixture, character- 
ized by a Gaussian field statistics. Any correlation 
function of arbitrary order is thus expressed via the 
second-order correlation function 



r(x,x') = (anx)a(x')), 



(1) 



where a denotes the boson field operator of the 
speckle beam. Notice that we are dealing with 
classical fields, so that the field operator a could 
be replaced by a stochastic c-number field, and the 
quantum averages by statistical averages over in- 
dependent data acquisitions. However, we prefer 
to keep a quantum formalism in order to outline 
the parallelism with the quantum entangled beams 
from PDC. 

3. Information about the object is extracted by mea- 
suring the spatial correlation function of the inten- 
sities (/i(xi)/2(x2)), where 1 and 2 label the test 
and the reference beam, respectively, and /^(xi) are 



the operators associated to the number of photo 
counts over the CCD pixel located at x^ in the 
ith beam. All the information about the object 
is contained in the correlation function of intensity 
fluctuations, which is calculated by subtracting the 
background term (/i(xi))(/2(x2)) : 

G(X1,X2) - (/I(xi)/2(X2)) - (/I(xi))(/2(X2)) . (2) 

By using the input-output relations of the BS, and 
the standard properties of Gaussian beams, the 
main result obtained in Ref. Il3l was 



G(xi,X2) = |rt|' 



ydx;ydx^/lt(xi,x'i)/i2(x2,X^)r(x;,x'2) 



(3) 



where /ii, /12 are the impulse response function de- 
scribing the optical setups in the two arms, and r, t 
the reflection and transmission coefficients of the 
BS. 

Equation (Q has to be compared with the analogous 
result obtained for PDC beams 12] : 



Gpdc(xi,X2) 



/ dx'i / dX2/li(xi,xi)/l2(x2,X2)rpdc(xi,X2; 



(4) 



where 1 and 2 label the signal and idler down-converted 
fields fli, 02, and 



rpdc(xi,X2) = (ai(xi)a2(x2)) 



(5) 



is the second-order field correlation between the signal 
and idler (also called biphoton amplitude). 

Thus, ghost imaging with correlated thermal beams 
presents a d eep analogy with ghost imaging with entan- 
gled beams [i3i[l3;L27]: they are both coherent imaging 
systems, which is crucial for observing interference from 
a phase object, and they offer analogous performances 
provided that the the beams have similar spatial coher- 
ence properties. They differ in (a) the presence of ft.* at 
the place oi hi (which in our case turns out to have no 
implications) and (b) the visibility, defined as 



V 



G 



(hh) 



G 



{h){h 



G 



(6) 



In the thermal case G(xi,X2) < (/i(xi))(/2(x2)) so that 
the visibility is never above i. Conversely, in the PDC 
case it can be verified that Gpdc/{h){h) scales as l-|--7iy, 

re (n) 
Ref. O) 

(n) ^ 1, the visibility can be close to unity, while bright 
entangled beams with (n) ^ 1 show a similar visibility as 
the classical beams. However, despite never being above 
i in the classical case, we have shown [ij,[i3,|23] that the 
visibility is sufficient to efficiently retrieve information. 



where (n) is the mean photon number per mode (see, e.g., 
j). Only in the coincidence-count regime, where 



The result of a specific correlation measurement is ob- 
tained by inserting into Eq. © the propagators describ- 
ing the setup. In the case of the ghost diffraction scheme 

of Fig. m /ii(xi,x;) = (zAF)-ie--i-i'=/^robj(x;), 



where Tobj(x) is the object transmission function, k 



27r/A, and /i2(x2,x'2) ^ {iXF)-^e- 



<k/F^ We get 



G(xi,X2) (X y'dx;y'dxie<''--'i-''--^)'=/^T4j(xl)r,(xl,x^) 

2^ /"dx'.F{robj} [(Xi - X2 - X')fc/F] rf(x2,X2 + X'] 



(7) 
(8) 



r 



where r„ and Ff denote the second order field cor- 
relation function defined by Eq. ^, as measured at 
the object near-field and far-field planes, respectively; 
^{Tohj} (q) = / ff e-^i-'^Tobjlx) is the amplitude of the 
diffraction pattern from the object, and Eq. is ob- 
tained from Eq. with some simple passages. 

We notice first of all that the result of a correla- 
tion measurement is a convolution of the diffraction pat- 
tern amplitude with the far-field coherence function Ff. 
Hence the far- field coherence length (denoted by Aa^f) 
determines the spatial resolution in the ghost diffraction 
scheme: the smaller the far-field speckles, the better re- 
solved is the pattern. In the limit of speckles smaller 
than the scale of variation of the diffraction pattern, 
we can approximate the far field coherence function as 
Ff(x2,X2 -f x') « (5(x')(/f(x2)), where (/f(x2)) is the in- 
tensity profile of the input speckle beam as observed in 
the far field [notice that in our ghost-diffraction setup 
(/f (X2)) = (/2(x2)), a part for some trivial proportional- 
ity factor]. In this limit we get 

G(xi,X2) ex |.F{Tobj}[(xi-X2)fc/F]|' af(x2))', (9) 

which means that the diffraction pattern of the object 
can be observed in the correlation function, when this is 
evaluated as a function of X2, for fixed Xi. The diffraction 
pattern is modulated by (/f (x2))^: hence, in order to ob- 
tain the whole diffraction pattern, the far-field intensity 
distribution must be sufficiently broad, so that (/f(x2)) 
is nonzero in the region where the diffraction pattern is 
nonzero. It can be easily seen that this condition turns 
out to be equivalent to requiring spatial incoherence of 
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the speckle beam illuminating the object, that is, the 
near- field coherence length (denoted here as Axn) must 
be small as compared to the object scale of variation. 

Eq. (O shows that the diffraction pattern can be also 
obtained by plotting the correlation as a function of the 
distance between the two points. By fixing this distance 
as r = Xi — X2, and varying the pixel positions in both 
arms as Xi and X2 = Xi — r we perform a spatial average 
of the correlation function, which amounts to measuring 

MM 

/"dxG(x,x- r) ex |.F{robj} [rfc/Fll^y dx(/f(x- r))2 . 

(10) 
If the spatial average is performed over large enough re- 
gions, the integral on the right hand side does not de- 
pend on r and is a constant. As already pointed out in 
Ref. [3, in this case there is no need of demodulating 
the correlation by the mean intensity in order to obtain 
the diffraction pattern. 

Second of all, and most important for the results pre- 
sented here: since the Fourier transform of the amplitude 
of the object transmission function is involved in Eq. (|SJ), 
the imaging scheme is coherent despite the fact that the 
beams are incoherent. Thus, ghost diffraction of a pure 
phase object can be realized with spatially incoherent 
pseudo-thermal beams. 

Quite different are the results for the HBT-type 
scheme, where the BS is effectively placed after the 
object. In this case the reference kernel changes to 
/^2(x2,x'2) = /ii(x2,x'2) = (zAF)-ie-— i'^/^Tobj (x^) , 
and a result of correlation measurement gives 



Ghbt(xi,X2) cx 



|dx;|dx'2e^^-/^(--'^---0r:,j(x;)robj(x^)Fn(xi,xi) 



In the limit of spatially incoherent light Fn(x'j^,X2) -^ S{x[ — X2)(/„(x']^)), and Eq. Hll|l can be re-casted as: 



Ghbt(xi,X2) cx 



/"dx'.F{|robjP} [(X2 -Xi +x')fc/F]Ff(x2,X2 +X') 



(11) 



(12) 



By comparing with the result of Eq. © for ghost diffrac- 
tion, we see that the HBT scheme only gives information 
about the Fourier transform of the modulus squared of 
the object transmission function: in the limit of incoher- 
ent light the imaging scheme is incoherent. Thus, the 
phase information about the object is lost and the HBT 
scheme is able to see interference only from absorptive 
objects. 

We can now consider the opposite limit, of spatially 
coherent light illuminating the object, achieved when the 
coherence length Axn (the speckle size) at the object 
plane is large compared to the object size. In this case, 
the HBT scheme allows to retrieve the diffraction pattern 
even of a pure phase object |35j- In this limit, in fact, the 
coherence function r,j can be taken as roughly constant 
over the regions of integration in Eq. Hll|) . which hence 
gives: 



Ghbt(xi,X2) (X |.F{Tobj}(xifc/F)|2 
x|.F{robj}(x2fc/F)p. 



(13) 



Evidently, if we fix xi and evaluate the correlation as a 
function of X2 we observe the diffraction pattern, even of 
a pure phase object. 

Notably, in the ghost diffraction scheme no diffraction 
pattern appears in the correlation as a function of the 
reference pixel position X2 for spatially coherent light. 
In the limit of spatial coherence, Eq. {TJ factorizes into 
the product of two integrals, one showing the diffraction 
pattern of the object in arm 1, as a function of xi, the 
other showing the mean far-field intensity profile in arm 
2. This can be readily seen also from Eq. ||SJ), where the 
limit of spatially coherent light at the object plane (limit 
of a single large speckle illuminating the object) amounts 

to (/f(x2)) ^ (5(X2). 

In practice, the classic HBT scheme uses the cross- 
correlation of two beams split on a BS after the object 
as a convenient way of measuring the auto-correlation of 
the beam transmitted through the object, as, e.g., done 
in Ref. J7| . We will actually measure the auto-correlation 
of the light in the test arm in the focal plane of the lens, 
defined as 

Cauto(x,x') = (/l(x)/i(x')) - (/i(x))(/i(x')) . (14) 

A part from a small shot-noise contribution at x' = x 
and some irrelevant proportionality factors, the results 
expected from such a measurement coincide with those 
of the HBT schemes described by Eq. (|11|) . and in the 
proper limits by Eqs. (fT^ . ((T^ . 

This comparison between the HBT and the ghost- 
diffraction schemes indicates that the measurement of the 
autocorrelation function is a valid test to prove whether 
or not a given object alters the phase of the incoming 
light or not: in the presence of spatially incoherent light, 
the ghost diffraction scheme fully preserves the informa- 
tion about the object phase in the diffraction pattern, 
whereas this information is lost in the autocorrelation. 
For a pure phase object, no interference pattern at all 



should appear in the autocorrelation. However, as it will 
become clearer in the next sections, the autocorrelation 
is extremely sensitive to the degree of spatial incoherence 
in the beam: even a small partial coherence in the incom- 
ing beam is sufficient to preserve some phase information 
in the autocorrelation function. 

We finally stress that despite having discussed the 
ghost diffraction results in the framework of classical 
speckle light [i.e., for which Eq. Q holds], the results ((HJ- 
(|13ll hold also for the entangled case [for which Eq. (0 
holds]. 



B. The object 

We have chosen a commercially available TGBS as our 
pure phase object. It is well known that such a device 
transmits the incoming light (and has close to zero ab- 
sorption/reflection) so that in the far zone of the object 
several distinct peaks are observed. This is because the 
diffraction angle obeys the thin grating equation 



i) = n\/d, 



0,±1,±2,.. 



(15) 



This equation holds when the incoming light is a plane 
wave normally incident on the grating, and tells that the 
light is observed in several orders n of the diffraction an- 
gle 9m and that the location of the diffraction peaks is 
found according to the ratio \/d where A is the wave- 
length of the light and d is the period of the grating. 

When observing the intensity distribution of the trans- 
mitted light in the far zone, the strength of the nth 
diffraction peak is rjn = |c„|^, where c„ are the diffrac- 
tion coefficients. Since a grating can be thought of as an 
infinite repetition of a single period of the grating Tobj , 
we may use Fourier series theory to write the diffiaction 
coefficients )32j 



1 



dxTobj(x)e-"2--/rf. 



(16) 



A TGBS is a grating made of a completely transpar- 
ent material which has grooves cut on the exit side so 
that a phase difference is imposed on the field exiting 
the grating depending on the position. This phase dif- 
ference can be expressed through the groove depth 5 as 
A0 — {ug — I)27r5/A, where Ug is the refractive index of 
the grating material. Thus, 5 can be chosen to give the 
desired phase shift. 

Typically, square gratings of width a (and period d) are 
used for TGBS, so that within the period d the object 
transmission function is 

^"•^^(^^^I 1, a<x<d- ^^^^ 

Calculating the diffraction coefficients gives 

% = l + sin2(A(/./2)^(a/d-l) (18) 

d 

f]n = sin^(A(j)/2)—-sinc^(nna/d), n = ±l,±2,(19) 
d^ 
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2 

x^ (mm) 

FIG. 2: Numerical simulations demonstrating the transi- 
tion between coherent and incoherent illumination of the ob- 
ject. Average intensity (Ji) in the far field of the test arm 
for different values of the speckle size at the object plane 
Ax-n = 2, 6, 12, and 48 /xm Numerical parameters: we use 1024 
grid points, and 32 pixels per groove period; d = 12.5 pm; 
a — 4.2 fim, Atf) — 0.847r; Ax{ — 80 ^m; averages are done 
over 10"" independent realizations. 



Choosing A(j) and the ratio a/d properly one can engineer 
the peaks to have the desired distribution. 

Our TGBS is from Edmund Optics and has 80 grooves 
per mm (stock number NT46-069), i.e., d = 12.5 /im. It 
is designed for A — 633 nm to have 25 % of the power in 
the and ±1 peaks, and 5 % in the ±2 order peaks. This 
means rjo — r\\ and rii/rj2 — 5, implying that At/) = O.TItt 
and that a/d = arccos(l/-\/5)/7r giving a — d/2.84 — 
4.4 /im. The smallness of a sets a limit for the relative 
coherence of the object beam, as discussed in detail in 
what follows. 

Since we are using a frequency doubled Nd:YAG laser 
with A = 532 nm the phase difference is instead A0 = 
flfO.TlTT ~ 0.847r. This implies 771/770 = 2.2, i.e., the 
central peak should be roughly a factor of 2 weaker than 
the ±1 order peaks when using light at this wavelength. 



C. Preliminary investigation through numerics 

In this section we will discuss some numerical results 
we used to select the size of the speckles we would need 
in the object plane in order to render the beam spatially 
incoherent. The numerical simulations were done using 
the method described in Ref. |22,] . Essentially, we Fourier 
transform noise convoluted with a Gaussian to obtain a 
certain speckle size in the object plane Axn (the size of 
which is controlled by the waist of the Gaussian). Impos- 
ing a pinhole of size iPph on this field and performing an- 
other Fourier transform gives a speckle field as observed 
in the far field of the object plane; the speckle size there 
Axf is determined by I^ph and is therefore controlled 
independently of Axn. The object used in the simula- 



tions was a purely transmissive square grating chosen to 
mimic the object predicted theoretically in Sec. lIIIi^ see 
Fig. 13 for specific parameters. Since ghost imaging im- 
plies that coherent imaging is done with incoherent light, 
it is important that the light is really incoherent relative 
to the object details. Thus, in a direct observation of 
the test arm far field average intensity, (/i), we should 
not be able to see the diffraction pattern of the object. 
As shown in Fig. 13 when Axn ~ 2^m (/i) does not re- 
veal any information about the diffraction pattern. Thus, 
such a speckle size corresponds to practically incoherent 
illumination of the object. However, as the speckle size 
is increased (Axn = 12 /im) (/i) reveals more and more 
information about the diffraction pattern, correspond- 
ing to partially coherent illumination. For large speckle 
sizes (Axn = 48 /im) (/i) is very close to the analytical 
diffraction pattern and the illumination is close to being 
completely coherent J36j . 



As shown in Sec. IIII Al if the test beam is incoher- 
ent with respect to the object, the autocorrelation of /i 
should not reveal at all the diffraction pattern of a pure 
phase object, while the pattern should appear in the 
cross-correlation. Fig. |3 compares the cross-correlation 
function G'(xi,X2) and the autocorrelation Cauto (xi , x'j^ ) 
of the test arm intensity, for different sizes of the near- 
field speckles. Both functions are calculated for a fixed 
xi. For Aa;,i = 2.0 /xm the cross-correlation shows a 
diffraction pattern that, as we verified, is very close to 
the pattern analytically calculated. On the contrary, the 
autocorrelation shows very little information in the side- 
bands; the n — \ sideband is 2.4% of the central {n = 0) 
peak. This confirms that this speckle size corresponds to 
practically incoherent illumination. The fact that some 
information is observed in the sidebands at all is because 
the 2.0 /im speckle on the scale of the object is not vanish- 
ingly small, but merely small. By increasing the coher- 
ence of the hght, we notice that the object information 
disappears from the cross-correlation and appears in the 
autocorrelation, in agreement with Eq. H13|l and with the 
discussion in Sec. IIII Al For Azn = 12 /xm the sidebands 
present in the autocorrelation coincide almost completely 
with the analytical diffraction pattern, showing thus that 
the autocorrelation function is sensitive to the presence 
of even a small partial coherence of the light. 



To conclude this discussion, we have to choose a 
speckle size at the object plane around 2 /zm in order 
to have beams that are truly incoherent with respect to 
the object details, so that information about the object 
is revealed neither in the autocorrelation nor the average 
of the test beam far-field intensity. 
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FIG. 3: Numerical simulations demonstrating the transition between coherent and incoherent illumination of the object, (a) 
The cross-correlation (solid line) normalized by the square of the mean intensity in the reference arm (the latter is shown by 
the dotted line) as obtained for different sizes Aa^n of the near field speckles, (b) The corresponding autocorrelation functions 
of the test arm intensity. The black squares indicate the peak values of the analytical diffraction pattern. Parameters: as in 
Fig.H 



IV. EXPERIMENTAL RESULTS 

A. Pseudo-thermal source and speckles from 
near-field scattering 

As shown in the previous section, the physical size of 
the TGBS requires a speckle size at the object plane Axn 
less than the finest object detail, i.e., less than 4.4 /im. 
Such a small size is not easy to achieve, but we managed 
to create speckles with Axn = 2.0 /im by placing the ob- 
ject very close to our pseudo-thermal source (z = 18 mm, 
practically limited only by the physical size of the half- 
inch cube BS). The speckles generated this way are the 
so-called near-field speckles (NFS) [Sfl, |31|j which are 
remarkably different from the classical far-field speckles 
(FFS), whose size is determined by the well-known Van 
Cittert-Zernike theorem |33 |. 

This part of the setup is quite different from what we 
used in our previous experiments in Refs. |17l llfll l27l | 
where the object was in the far zone of the source. In- 
stead, in the current setup the source is very close to the 
object, and, as explained in detail below, at a given point 
of the object plane the waves interfering do not originate 
from the entire illuminated spot. 

To understand how NFS are formed, let us first de- 
scribe the pseudo-thermal source in some detail. 

Our thermal source consists of a laser beam illuminat- 
ing a slowly rotating ground glass, followed by a square 
cell 5 mm thick, containing a concentrated solution of la- 
tex particles with average diameter p = 3 /im. The cell is 
almost in physical contact with the rotating glass and on 



its exit face there is a pinhole (diameter Dp\^ ~ 4.5 mm) 
which determines the transverse size of the source. The 
combination of the ground glass and the turbid solution 
is an easy and convenient way to generate truly ran- 
dom speckles. Indeed, the ground glass alone would pro- 
duce only partially stochastic speckles because the pat- 
tern would be reproduced after one turn of the glass. On 
the other hand, the turbid cell guarantees stochasticity, 
but if used alone would exhibit a residual transmitted 
coherent component which is clearly undesired. 

Due mainly to Brownian particle motion (and secon- 
darily to glass rotation), the speckle pattern generated 
by the source fluctuates randomly with time and is char- 
acterized by a coherence time Tcoh which can be tuned 
by varying the turbidity of the solution. In our case we 
had Tcoh ~ 10 ms. 

Multiple scattering occurs inside the cell, so that the 
light beam exiting the source has a divergence angle 8eff 
larger than the angle that would be expected from single 
scattering. The latter is associated to particle diffrac- 
tion and given by O^if ~ ^/ P- The effective value of 
Ocff depends on the detailed features of the scattering 
cell, as particle concentration and cell length; in prac- 
tice, we can claim that our thick thermal source behaves 
as an ideal thin thermal source characterized by spa- 
tial inhomogeneities or " scatterers" of effective diameter 

PcS « P (Bdif/Oeff) < p. 

When the light generated by such a source is observed 
at a plane located at a sufficiently large distance z, each 
point of this plane is reached by the contributions emerg- 
ing from the entire radiating region ZJph. Under this 



condition, the stochastic interference between the (many) 
different waves gives rise to a speckle- hke intensity distri- 
bution {far-field speckle, FFS), whose correlation function 
is described by the Van Cittert-Zernike theorem, and is 
characterized by the average speckle size |33|| 
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Thus, the requirement for obtaining FFS is QetfZ ^ i'ph, 
i.e. 
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When this condition is not fulfilled, the waves interfering 
at each point of the observation plane originate from a 
region D* « 0eff z « ^- z smaller than the radiating 
region -Dph. Provided that D* is not too small, at a given 
point in the observation plane one would still get contri- 
butions from many different scatterers, which is sufficient 
to produce near-field speckles. Applying again the Van- 
Cittert Zernike theorem with D* instead of Z3ph, one gets 
the remarkable result Isfl 



Ax « peff (NFS) . 



(22) 



according to which the average NFS size is only deter- 
mined by the effective size of the scatterers, and is inde- 
pendent of both A and z. To fulfill the criteria for NFS 
we must have (i) D* < D and (ii) many scatterers inside 
D*, e.g., {D* / pcsY > 1- This imphes that the distance 
z must fulfill the two conditions 
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In our experiment we have A = 0.532 //m, Z?ph — 4.5 mm 
and pcff < 3 /ini. Thus the second criterium is easily 

2 2 

fulfilled because ^^ < ^ = 17 /im, so that z 3> 17 /im 
is enough. The first criterium is somehow more difficult 
to evaluate, because it requires more detailed knowledge 
of pcff • However, our final purpose was to make speckles 
as small as possible and, therefore, we set the distance 
between the object and the pseudo-thermal source to its 
minimum value z = 18 mm 37|, which to a large extent 
meets the first criterium in Eq. (|23|l . 

Then we measured the speckle pattern in the object 
plane by removing the lens F and inserting a x20 objec- 
tive to image the object plane on the CCD (the magnifi- 
cation is needed because the CCD pixel size is 6.7 ^m). 
The speckle size was finally estimated by performing the 
autocorrelation of such pattern as shown in Fig. ^ The 
peak above the baseline had a full-width half maximum 
(FWHM) of 1.98 ± 0.02 /^m. This gives an estimate of 
the near-field speckle size (see Refs. [TJjIlJ for more de- 
tails) Aa;,! « 1.98 /im, which, as shown in Sec. IIIICI 
should be small enough to render the beams incoherent 
with respect to the object details. 

For completeness, we also measured the size of the 
speckles in the far-field (the size of the speckles on the 
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FIG. 4: Spatial autocorrelation function of the speckle beam 
just before the object (/n(x')Jn(x' -I- x)). The peak above 
the baseline has a FWHM of 1.98 ± 0.02 /im, which gives the 
characteristic size of the near field speckles. 



CCD) by removing the objective and reinserting the lens 
F. The procedure gave Axf — 11.1 ± 0.1 /im. This value 
actually overestimates the real size of the far-field speck- 
les, because the CCD pixel size (6.7 /im) is too large with 
respect to the speckles, so that the speckle pattern under- 
goes a substantial smoothing. In any case, the measured 
Aa;f determines the spatial resolution of the ghost diffrac- 
tion pattern, and in our case turns out to determine the 
width of the ghost diffraction peaks. 



B. Ghost diffraction versus HBT scheme: case of 
incoherent illumination 

We performed a first set of measurements keeping the 
object plane at a distance z = 18 mm, thus having speck- 
les at the object plane of size Ax^ ~ 2 /im. 

In order to characterize the diffraction pattern created 
by the TGBS and to provide a reference for the ghost 
diffraction pattern, we performed preliminary measure- 
ments with coherent laser light: we removed the scat- 
tering media from the setup of Fig. ^ and recorded the 
transmitted light of the TGBS in the focal plane of the 
lens F . This measurement was already not straightfor- 
ward because of the large values of the scattering angles 
of the grating equation <fT5|l : the nth order peaks at an- 
gles On = nX/d, are displaced in the far field at positions 
Xn = 0nF = nFX/d. Using the numbers in our setup 
we have x±i = ±2.13 mm and x±2 = ±4.26 mm, so 
that the distance between the two 2nd order peaks is 
larger than the extension of our CCD. We therefore had 
to do 3 measurements in order to observe all the peaks: 
first n — —2,-1, and were observed, then the CCD 
was shifted to observe n = —1,0, and +1 and finally 
ri = 0, +1, and +2. A second problem was represented by 
the small width of the diffraction peaks (few pixels on the 
CCD) which provided a too poor sampling of the diffrac- 
tion pattern. In order to evaluate the relative heights 
of the peaks we performed an integration in the region 
around each peak, which gave r]i/r]o = 4 and rji/rj2 = 2. 
Moreover the diffraction pattern is not symmetric and for 
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FIG. 5: Experimental demonstration of ghost diffraction of a 
pure phase object using incoherent classical light, (a) Snap- 
shot of the speckles recorded by the CCD in the far-field 
plane. The reference beam is in the upper region and the 
white frame shows the region used for the correlation. The 
test arm is in the lower region, and the white symbols indi- 
cate the 3 different single-pixel positions used for the corre- 
lation, (b) Ghost diffraction patterns reconstructed via the 
cross- correlation between the test and the reference arm, mea- 
sured by locating xi at each of the 3 pixel positions, and by 
varying X2 (18000 averages). The 2-D plots of the cross corre- 
lation G(xi,X2)/(/f(x2))^ are shown as functions of X2 — xi, 
together with their 1-D cut in the horizontal direction. 



example 77-1 =/= rji. This is somewhat different from what 
we expected from the theory, and probably depends on 
some defects of fabrication of the TGBS, but it will serve 
as our reference. Experimentally we found the peaks to 
be located at x±i = ±2.22 mm and x±2 — ±4.43 mm, in 
good agreement with the theory. 

We also used the coherent illumination to set the ori- 
gin of our coordinate systems: in the test arm Xi = 
corresponds to the the n = diffraction peak, while in 
the reference arm the X2 = point is the location of the 
reference beam. If the test-arm pixel is for example fixed 
at Xi = in the subsequent correlation measurements, 
we expect that the ghost diffraction pattern will emerge 
in a region of the reference arm centered around X2 = 0, 
while by shifting Xi the pattern will shift accordingly. 

We then reinserted the scattering media to measure 
the ghost diffraction pattern of the TGBS from the cross- 
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FIG. 6: Quantitative comparison between ghost diffraction 
and laser illumination. The height of the diffraction peaks of 
the cross correlation shown in Fig. |3 are evaluated by per- 
forming an integral over the peaks. They are compared with 
peak values measured via coherent laser illumination. 



correlation between the two arms. A typical snapshot of 
what we observe on the CCD is shown in Fig.|Sfa). The 
upper part contains the reference arm intensity. For the 
correlation, we selected a narrow strip (128 pixel wide) 
centered around y — and extending over the entire 
X axis (1024 pixels). In the test arm no spatial infor- 
mation was extracted since there we collected the light 
from a single fixed pixel. Initially we located the pixel 
in the test arm at the point xi = [pixel at position 
P2 in Fig. ISJa)], and we measured the cross correlation 
(/i(xi = P2)l2{^2)) as a function of X2 varying in the 
region shown by the white frame in Fig. |5ja) , by aver- 
aging over 18000 snapshots. As seen from Eq. lO this 
gives a diffraction pattern that is centered on X2 = 0. In 
this case we are able to reconstruct only the n = —1,0 
and +1 peaks, because the higher order peaks are out- 
side of the reference region imaged on the CCD. In order 
to reconstruct also the n — ±2 peaks, we repeated twice 
the measurement by shifting the test arm pixel at the 
positions Pi and P3, respectively, [see Fig. [Sfa)], so that 
the the diffraction pattern emerging from the correlation 
shifts accordingly, as dictated by Eq. ©. 

The result of these measurements are shown in 
Fig. Elb), which plots the cross correlation scaled to 
(/f(x2))^. The 2D reconstructed diffraction patterns are 
shown close to their cut along the x-direction, for each 
positioning of Xi . Since each diffraction peak covers only 
few pixels, the x-cuts of Fig. |SIb) only give a qualita- 
tive image of the diffraction pattern, but do not allow a 
quantitative estimation of the peak heights, due to the 
poor sampling. In order to extract the relative height 
of the peaks, we located groups of pixels having a sub- 
stantial value above the noise floor and added together 
their values, which effectively corresponds to integrating 
over each peak. The results are shown in Fig. |S1 and 
compared with those obtained from coherent laser illu- 
mination with the same technique. We observe that they 
agree extremely well: we have successfully created the 
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FIG. 7: Cross-correlation function (/2(x2)/i(xi)), calculated 
with the spatial average technique. The upper part is a 2D 
plot of the correlation function (100 averages). The lower 
plot displays the integrals over the peaks compared with the 
results of laser illumination. 



correct diffraction pattern of the pure phase object from 
the correlation. Also notice that the n = ±1 peaks are 
recorded twice, and the n = is recorded three times: 
these overlaps happen when displacing the single pixel 
position, and they agree very well with each other. Fi- 
nally, note that the position of the diffraction peaks agree 
well with the theoretical prediction in Sec. IIIIBI 

As discussed in Sec. lIII Al the spatial average technique 
provides a a faster and more efficient way of measur- 
ing the ghost diffraction pattern. In this case the cross- 
correlation is measured by varying both X2 and Xi for a 
fixed X2 — xi . The results of this kind of measurements 
are shown in Fig. |7| where the upper part of the fig- 
ure is a 2-D plot of the measured cross-correlation, and 
the lower part displays the integral over the diffraction 
peaks compared to the laser illumination results. Also 
in this case the agreement is excellent. From this fig- 
ure we notice that only few averages over snapshots are 
needed to reconstruct the diffraction pattern, because a 
large number of averages over spatial points is performed, 
thus increasing the convergence rate [1^, H^ ■ Moreover 
the whole diffraction pattern is reconstructed in a single 
measurement. Despite being much more efficient than 
the single-pixel reconstruction, the spatial average tech- 
nique does not follow the ghost imaging original spirit, 
which assumes that the imaging information is extracted 
by only operating on the reference arm. In this case, spa- 
tial information is also extracted from the test beam 1, 
by varying the pixel xi . 

As a straightforward demonstration of the degree of 
incoherence of the beams used, we present in Fig. |H1 a 
measure of the autocorrelation function in the test arm 
C'auto(xi,x']^). This is measured by fixing x.[ at position 
P2 and by varying Xi, in the same way as described in 
Sec. nil CI Evidently it does not reveal any significant 
information about the diffraction pattern. In fact, the 
first order peaks are barely visible and are at a level of 
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FIG. 8: Autocorrelation function (Ji(xi)/i(xi)) of the test 
arm, measured by fixing xi at position P2 and varying xi 
(10000 averages). 



8% of the main peak, in trend with the prediction of the 
numerical results of Sec. IIII(]I 



As argued in Sec. IIII Al 
this type of measurement is equivalent to a HBT-type 
scheme, and it works as an incoherent imaging scheme 
when using incoherent light; in this case it is expected 
to give no information about a pure phase object. We 
can thus conclude that i) the TGBS is truly a pure phase 
object and ii) the speckle light we use is truly incoherent 
relative to the object. 



C. Ghost diffraction versus HBT schemes: case of 
partially coherent illumination 

In this section we present results obtained by gradually 
increasing the spatial coherence of the light illuminating 
the object. This is achieved by increasing the distance z 
between the pseudo-thermal source and the object plane 
(see Fig. Pi . 

We performed a second set of measurements with this 
distance set as z = 115 mm. The measured autocorrela- 
tion of the light illuminating the object gave a speckle size 
Axn — 14 fj,in (FWHM of the autocorrelation peak). The 
main results obtained in these conditions are displayed in 
Fig. El In a third set of measurements the object-source 
distance was z = 300 mm, and the measured speckle size 
was Axn = 33 fiva. Figure [TUI displavs the results in this 
case. 

By increasing the source-object distance the light gains 
some partial coherence relative to the object. This is al- 
ready evident by the distribution of the speckles recorded 
in the far field of the test arm, shown in the lower parts of 
frames (a) in Figs. 1^1 and 1 101 Differently from the case of 
incoherent illumination, where the mean intensity distri- 
bution of the test arm is almost flat [see Fig. [Sfa)] , two 
broad peaks in correspondence of the n — ±1 diffrac- 
tion orders are now clearly distinguishable in the speckle 
distribution. As the coherence of the light increases 
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FIG. 9: Case of partially coherent illumination. The distance 
source-object is z = 115 mm, with Axu = 14 /im. (a) Far- 
field speckle distribution in a single snapshot (upper part: 
reference arm, lower part:test arm), (b) Cross correlation of 
the test and reference arm, for xi fixed at the origin and 
30000 averages. The dashed curve shows the reference in- 
tensity squared, (c) Autocorrelation of the test arm light 
intensity. Full line: horizontal section of the autocorrelation 
function. Open triangles: peak values of the autocorrelation 
(integral over the peaks). Circles: peak values measured with 
coherent laser illumination. 



[Fig. llOf a)] they become narrower and more pronounced. 
Notice that the order peak is barely visible in these 
plots because its intensity is lower (as dictated by the 
TGBS). 

The cross-correlation between the test and the refer- 
ence arm, obtained by fixing Xi at position P2 , is plotted 
in frames (b) of Figs.l^andllOI. We see that by increasing 
the coherence of the light, the height of the ±1 diffraction 
peaks decreases with respect the order peak, and the 
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FIG. 10: The spatial coherence of the light illuminating the 
object is further increased: z = 300 mm, Axn = 33 /im. 
Frames (a)-(c) display the same quantities as in Fig. |n| 



diffraction pattern gradually disappears from the cross- 
correlation. Notice that these plots shows the "bare" 
cross-correlation function G(xi,X2) [i.e., there is no scal- 
ing factor (/f(x2))^]. As predicted by Eq. (Q, the corre- 
lation scales with the square of the mean intensity of the 
reference arm, whose profile is plotted by the dashed lines 
in the figures. By increasing the near-field coherence, the 
far-field intensity spot becomes narrower, until the mean 
intensity vanishes in the region where the higher order 
diffraction peaks should emerge. In principle, the cor- 
rect height of the diffraction peaks could be recovered by 
dividing the correlation by (/f(x2))^, but this operation 
also amplifies the noise in the regions where the intensity 
level is low. This is evident when the cross-correlation is 
normalized, as shown in Fig. ^] (see also Fig.O, and we 
notice that in the case of Axn = 14 /xm the ±1 peaks can 
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FIG. 11: Cross-correlation function normalized to the square 
of the mean intensity in the reference arm. (a) z — 115 mm, 
Axn ~ 14 /im: the ±1 peaks can still be reconstructed, (b) 
z — 300 mm, Axn = 33 fim: the diffraction peaks disappear 
in the noise. 30000 averages are performed for both frames. 



be almost reconstructed, while for Axn = 33 /im they 
disappear in the noise. 

In other words, by increasing the coherence, the signal- 
to-noise ratio for the reconstruction of the higher order 
peaks decreases, and the information about the diffrac- 
tion pattern becomes less and less accessible. To make 
this argument more formal, we remind that the signal- 
to-noise ratio is proportional to the visibility defined by 
Eq. ©, as derived in Ref. ^. By using Eq. ((SJ), 
we can readily conclude that for G small the visibil- 
ity of the nth order peak, located at X2 — X2, is 
V « G(xi,X2)/[(/i(xi))(/2(x2))] « r?„(/2(x2))/(/i(xi)). 
Since the point xi is fixed at P2 where the test intensity is 
nonzero, the visibility, and hence the signal-to-noise ratio 
is proportional to the intensity in the reference arm. 

Conversely, by increasing the near-field coherence, the 
diffraction pattern gradually appears in the autocorrela- 
tion function of the test arm, displayed by frames (c) of 
Figs. El and ^| In these plots the diffraction peak val- 
ues measured via the autocorrelation (triangles) are com- 
pared to the values measured by coherent laser illumina- 
tion (circles). For Axn — 33 /im the partial coherence of 
the light is already enough to permit an almost perfect 
pattern reconstruction in the autocorrelation function. 

The results presented in this section evidence a clear 
complementarity between the ghost diffraction scheme 
and the HBT scheme, that will be further discussed in 
the next section. 

A final remark is the following: had we used the spa- 
tial average technique, some information on the diffrac- 



tion pattern would have been preserved in the correlation 
when increasing the spatial coherence. In this technique, 
in fact, the pixel position Xi in the test arm is scanned to- 
gether with X2 ; in this way, if some information is present 
in the test arm intensity distribution, this is retrieved 
from the correlation. By increasing the spatial coher- 
ence, the diffraction pattern becomes visible in the in- 
tensity profile of the test arm as shown by Figs. EJa), 
llOf a). and becomes also visible in the correlation as a 
function of Xi. But, obviously, as the diffraction pattern 
appears in the test arm, is not possible any more to speak 
about "ghost diffraction". 



V. CONCLUSION 

We have shown that coherent imaging with incoherent 
classical thermal light is able to produce the interference 
pattern of a pure phase object. This provides the ulti- 
mate demonstration that entanglement is not needed to 
do coherent imaging with incoherent light, not even in the 
case of a pure phase object. As our g rou p has pointed out 
in previous pubhcations [UlllllllllSIIlllilll, 
the only evident advantage of using entangled light might 
be that of obtaining a better visibility. 

A remarkable aspect of the present experiment is 
the degree of incoherence of the pseudo-thermal speckle 
beams used. In order to render the beams incoherent 
with respect to the object (a standard transmission grat- 
ing beam splitter with 80 grooves per mm), we had to 
create speckles which in the object plane had a size of 
2.0 ^m. This was made possible by exploiting the so- 
called near-field scattering [Sfl |3ll , in which the speckles 
are created so close to the source that their size is gov- 
erned solely by the roughness of the scattering medium. 

In such conditions of spatial incoherence, we have 
shown that no information on the phase object is present 
in the light outcoming from the object: neither the far 
field intensity distribution of the test arm nor its auto- 
correlation function (HBT scheme) reveal the diffraction 
pattern. This information is instead present in the cross 
correlation between the test arm and a reference arm that 
never passed through the object (ghost diffraction). Our 
results indeed evidence that, when trying to extract infor- 
mation on a pure phase object, there exist a clear comple- 
mentarity between the ghost diffraction scheme and the 
HBT scheme. In the HBT scheme the presence of a cer- 
tain degree of spatial coherence is the essential ingredient 
that permits to extract some phase information, and the 
information becomes more correct as the coherence in- 
creases. Conversely, the ghost diffraction scheme works 
as a coherent imaging scheme only thanks to the spatial 
incoherence of the light, and the more the light is incoher- 
ent, the better the information is reconstructed. These 
results contradicts what was indicated in the introduc- 
tion of Ref. [3, where the possibility of doing coherent 
imaging in a ghost imaging scheme employing splitted 
thermal light was ascribed to the presence of spatial co- 
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